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We study quantum antiferromagnets on two-dimensional bipartite lattices. We focus on local 
variations in the properties of the ordered phase which arise due to the presence of inequivalent sites 
or bonds in the lattice structure, using linear spin wave theory and quantum Monte Carlo methods. 
Our primary finding is that sites with a high coordination tend to have a low ordered moment, at 
odds with the simple intuition of high coordination favoring more robust Neel ordering. The lattices 
considered are the dice lattice, which is dual to the kagome, the CaVO lattice, an Archimedean 
lattice with two inequivalent bonds, and the crown lattice, a tiling of squares and rhombi with a 
greater variety of local environments. We present results for the onsite magnetizations and local 
bond expectation values for the spin-1/2 Heisenberg model on these lattices, and discuss the exactly 
soluble model of a Heisenberg star, which provides a simple analytical framework for understanding 
our lattice studies. 

PACS numbers: PACS numbers: 75.10.Jm, 71.23.Ft, 71.27. +a 



I. INTRODUCTION 

The long-wavelength universal properties of antiferro- 
magnets on bipartite lattices are well established Q. In 
this paper, we ask to what extent local properties - the 
dependence of the ordered moment or bond strength on 
local coordination - display systematic behavior. To do 
this, we consider quantum Heisenberg antiferromagnets 
on various two-dimensional lattices. These lattices are bi- 
partite, so that classical antiferromagnets form collincar 
Neel states at zero temperature; for quantum spins, this 
ordering persists, albeit with reduction of the order pa- 
rameter due to quantum fluctuations. A review of results 
obtained for uniform Archimedean lattices (in which all 
sites - but not necessarily all bonds - are equivalent) is 
given in Ref. Q . Generally, the size of the order parame- 
ter depends on dimensionality - Heisenberg Neel order is 
absent in d = 1 - and coordination, z: it is bigger for the 
square (z — 4) than for the honeycomb (z — 3) lattice. 

Upon increasing the coordination z the exchange 
'mean-field' acting on a spin at a given lattice site grows. 
The larger this mean-field, the less effective quantum 
fluctuations are in reducing the staggered alignment of 
the spin at that site. For inhomogeneous lattices with 
sites of different coordination, this suggests that the lo- 
cal staggered magnetization should increase with sites 
coordination z. Analogously, one would expect average 
bond energies to decrease with z, as fluctuations trans- 
verse to the ordered moment become increasingly difficult 
to coordinate between a growing number of neighbors. 

However, recent investigations on a two-dimensional, 
quasiperiodic system, the octagonal tiling, with six dif- 



ferent local coordinations, z — 3,4,5,6,7,8 lead to the 
result that the local staggered magnetization does not 
follow this expectation; instead, both the local staggered 
magnetization rtii and the averaged local bond strength 
of all bonds connected to a given site, tend to de- 
crease with increasing z (cf. Sec. II for a formal defi- 
nition of both rrii and hi). This property has been ob- 
served in quantum Monte Carlo simulations Q , in linear 
spin wave theory £| and in a real space renormalization 
group approach based on the self-similar structure of the 
quasiperiod lattice The absence of translation in- 
variance, however, has important consequences for the 
spectrum and eigenmodes of the spin Hamiltonian Q, 
and it is the invariance under scale transformation that 
determines the real space properties of the antiferromag- 
net 0, 0| . In order to assess the relevance of the absence 
of translation symmetry on the behavior of the local or- 
der parameter distribution, here we consider various in- 
homogeneous but periodic lattices. 

In particular, we calculate the variations in mi and b t 
for lattices with small unit cells that contain inequivalent 
bonds and/or sites. We consider the dice lattice, the 
CaVO lattice and a lattice that to our knowledge has not 
been studied previously, the crown lattice, which is built 
from squares and rhombi. This structure shows similar 
local environments to those found in the octagonal tiling, 
but unlike the latter possesses translational invariance. 

Our results on regular lattices are in qualitative agree- 
ment with those of the octagonal tiling. To provide a 
simple analytical explanation of the z-dependence of the 
local quantities, we appeal to the exactly soluble case 
of the Heisenberg star Q. In particular, this resolves 
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the conundrum of a decreasing rrii(z) by showing that, 
while the average moment ra™ increases with the average 
coordination z av of the lattice, an order-by-disorder type 
mechanism concentrates the reduction of the ordered mo- 
ment to highly coordinated sites. 



II. MODEL AND NUMERICAL METHODS 

In the following we consider the nearest neighbor 
Heisenberg antiferromagnet, described by the Hamilto- 
nian 

H = JJ2S 1 S J (1) 

(id) 



III. LOCAL PROPERTIES OF 
TWO-DIMENSIONAL BIPARTITE LATTICES 

III.A The Dice lattice 




where the sum is taken over all bonds in the lattice. J > 
is assumed to be independent of the local environments, 
and set equal to 1 in the rest of the paper. We employ 
a combination of linear spin wave theory Q and quan- 
tum Monte Carlo simulations in our analysis. Details of 
the linear spin wave approach will be presented in the 
following section for each specific lattice. 

The quantum Monte Carlo simulations wereperformed 
using the stochastic series expansion method |$| for finite 
lattices. In particular, for the dice lattice we considered 
systems with N s = 192, 432 and 768 sites, for the CaVO 
lattice those with 64, 144, 265 and 748 sites, and for the 
crown lattice, systems with 112, 252, 448 and 700 sites. 
In each case the temperature was taken sufficiently low to 
obtain ground state expectation values. For a given site 
i of the finite size lattice, the local value of the staggered 
magnetization is given by 
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FIG. 1: The dice lattice, showing a localized excitation. 



Our simplest example is the dice (or T3) lattice, Fig.^ 
which has a three-site unit cell for which Eq. (1) is easily 
diagonalized in linear spin wave theory. This is a system 
with a net magnetic moment per unit cell and therefore 
it has a Goldstone mode with a quadratic dispersion as 
in a ferromagnet. For completeness, we present an out- 
line of the spin wave calculation for this simple case. The 
magnetic and the structural cell unit cells are the same 
and have three sites, one on sublattice A with z = 6, 
and two on sublattice B, both with z — 3. One intro- 
duces Holstein-Primakoff operators of three types, cor- 
responding to each of these sites: ai,bi,Ci, (i = 1, ...,N 
where N is the number of unit cells) along with their ad- 
joints, obeying the appropriate bosonic commutation re- 
lations. Introducing Fourier-transformed operators — 



Hamiltonian is Hu 
N b = 2N, where 



i 

Vn ■ 
-JS(S 



lkVi b(c)i, the linearized 
- l)N b + JSHW, with 



where the sum extends over all the Ni lattice sites j which 
are equivalent to site i with respect to the antiferromag- 
netic unit-cell. Finally, standard finite size scaling analy- 
sis was performed to obtain the local staggered magneti- 
zation in the thermodynamic limit Q . For a given bond 
on the lattice, we define the corresponding bond 
strength as bij — \ (SiSj)\. Similarly, the averaged local 
bond strength at a site i is 



1 2 

bi = - ^2\{SiSj 



(3) 



where j extends over all the z sizes that are connected to 
site i |3j . In the following we present the results of the 
quantum Monte Carlo simulations as well as those from 
linear spin wave theory for the different lattices. 




(4) 



Here, A(k x ,k y ) = e ~ l ^ k -/ 2 {e ak y/ 2 + 2cos(-y/3fc a; /2)) 
(k x , k v are expressed in dimensionless units). This Hamil- 
tonian is diagonalized by a generalization of the standard 
case as described in Q in order to find a set of bosonic 
operators at, 7g such that 



^ (2) = >>^7 £ ) 



u>i{k) 
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Lu 2 (k) 

UJ 3 (k) 



.(5) 



An analytic solution is easily found and the results for 
the variations of u>i along different symmetry directions 
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are shown in Fig. |2J The flat branch 0*3 = 3 corresponds 
to an excitation that lives on a subset of sites, situated 
on hexagons of B-sites (as illustrated in Fig. In this 
linearized model, each of the six spins around the cen- 
tral site precesses in anti-phase with its neighbors. These 
states arise from the local topology, which gives rise, sim- 
ilarly, to the "Aharanov-Bohm cages" discussed in |9( for 
electrons subjected to a magnetic flux on the dice lat- 
tice. The ground state energy per site is E /N s —0.633 
within linear spin wave theory, compared to the quantum 
Monte Carlo result of -0.6384(1). 




FIG. 2: The three dice lattice energy bands along the main di- 
rections in k-space. Special points in the hexagonal Brillouin 
zone are O: origin, A: vertex, B: midpoint of an edge. 



The local staggered magnetization of A-sites, tua = 
S — (aj<Zj), and that of the B-sites, m^, can be calculated 
from the transformation matrices Og that take the set 
at, bt, ct into the set at, pt, 7!. One finds 



rriA 



TUB 



i - 1 E 

9 M 



1 - 1 E 

9 M £—1 



3 - ui 2 (k) 
•V^3 + 2w 2 (fc) 

3 — u>2(k) 
6 + Auj 2 {k) 



0.356, 



0.428. 



(6) 



The corresponding quantum Monte Carlo results are 
0.3754(3), and 0.4381(2), respectively. Note that the 
quantum correction to the classical B-sublattice magne- 
tization, Sm b = S — itib is half that of the A-sublattice 
correction 6m a = S — tua- The local magnetization on 
the dice lattice thus follows the trend seen on the oc- 
tagonal tiling, being smaller on sites with higher coor- 
dination number. As for the bond strengths, there is 
only one type of nearest neighbor bond in this system, 
b = \(SaSb)\ = 0.316 in linear spin wave theory (quan- 
tum Monte Carlo result: 0.3189(1)). 

III.B The CaVO lattice 

The second lattice considered here has four sites per 
unit cell, and twice that in the magnetic unit cell. This 



is the Til or CaVO lattice (after the calcium vana- 
dium oxide compound whose spins lie on the topolog- 
ical equivalent of this structure) discussed previously 
(c.f. discussion in Ref. Q). We introduce a number- 
ing from 1 through 4 for the four sites of a structural 
unit cell. The magnetic unit cell is doubled, with x 
and y axes rotated by 45° with respect to the original 
axes as shown in Fig. We thus introduce the boson 
destruction operators ai_4 (61-4) corresponding to sub- 
lattice A (B). The magnetic moment per magnetic unit 
cell is zero. The transformation using bosonic operators, 
linearization, and Fourier transformation for the sets of 
a,i and bi are carried out similarly as for the dice lat- 
tice. For a system of N magnetic unit cells, one finds 
Hn n = -JS{S + l)N b + JSHW, with N b = 12N and 



Hi 


H 2 




Hi 



(7) 



where H\ = zib{ 



35ij, H 2 is the adjoint of H 2 with 

4 o\ 



H 2 = 



1 
1 



(8) 
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where 



exp 



The basis vectors are 



(oi, a 2 , az,ai,b\,b\,b\, b\). Numerical diagonalization of 
this Hamiltonian leads to d^ = O^aj:, after discretiz- 
ing the Brillouin zone. We obtain four branches u>i 
(i = 1,..,4) shown in Fig. There are dispersionlcss 
directions, k x = ±k y for excitations living on strips, such 
as the one shown in Fig. [3 of energy E = 2V2. The par- 
ticipating sites occur with larger (smaller) amplitudes as 
indicated by filled (open) circles. 

The ground state energy per site for this system within 
linear spin wave theory, obtained by a discrete sum over 
the Brillouin zone, is Eq/N s m —0.5376 , in agreement 
with the value of Ueda et al. 0] . Ueda et al. also show 
that in this system, the Neel state corrected by spin waves 
has an energy very close to the energy of a plaquette res- 
onating valence bond (PRVB) state. The ground state 
energy obtained in linear spin wave theory is seen to com- 
pare less well with the quantum Monte Carlo results, 
—0.55367(2) than in the case of the dice lattice. Similar 
deviations between linear spin wave theory and quantum 
Monte Carlo results are also obtained for the local quan- 
tities: All sites in the CaVO lattice have the same value 
of the staggered magnetization, which is obtained within 
linear spin wave theory as to; — S — (ajaj) = S— (b\bi) = 
to. Our k-space sum extrapolates to m sa 0.213, close 
to the value, 0.212, published in [l^}, but significantly 
larger than the quantum Monte Carlo result, 0.1805(2). 
The increased relevance of quantum fluctuations, which 
are not captured appropriately within linear spin wave 
theory stems from the closeness of the system to a quan- 
tum critical point [llj , beyond which the system enters a 
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plaquette phase with dominant singlet formations on the 
squares in Fig. |3| This is reflected in the enhanced bond 
strength on the plaquettes, b p — 0.40731(6) compared 
to the remaining (dimer) bonds — 0.2926(1). Linear 
spin wave theory does not fully capture this effect, yield- 



b p = \ - S 2 



mg 
b d = \-S 2 



5(1 - 2m + 2Re(alb\)\ = 0.365, and 
- S(l-2m + 2Re(a\b\)\ = 0.344, respectively 




FIG. 3: The CaVO lattice, showing a state localized on a 
strip (open circles correspond to bigger amplitudes than the 
filled circles). 



En 




FIG. 4: The CaVO energy bands along the main symmetry 
directions. Special points in the square Brillouin zone are O: 
origin, X: midpoint of edge, F: vertex. 



III.C The crown lattice 

This is a square lattice system with seven sites per unit 
cell, labeled as shown in Fig. and twice this number in 
the antifcrromagnetic unit cell. Each unit has a mirror 
symmetry that reduces the number of distinct environ- 
ments to four. This periodic structure has most, though 
not all, of the local environments present in the quasiperi- 
odic octagonal tiling, and can in fact be obtained via pro- 
jection from a four dimensional cubic lattice using the 
method for obtaining approximants of the tiling outlined 
in [l^. Each of the sites occurs once on each sublattice, 
as for the CaVO lattice, and we introduce accordingly 
the Bogoliubov operators a, , aj for sites of sublattice A 




FIG. 5: A unit cell of the crown lattice, with sites numbered 
from 1 to 7. 



and bi, b\ for sites of sublattice B (i = 1, 7). The lin- 
earized Hamiltonian is H Hn = -JS(S + l)N b + JSH^ 2 \ 
where 2V& = 14N. H^> takes the form given in Eq. Q 
in a basis a T = (oi, ct2 ; &7i b\, b\, ■■, b\) where 



(-ffi)ij = Zi5ij\z= (6,3,3,3,4,4,5), 



(9) 
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(10) 



The 14 x 14 system is diagonalized to find the new ba- 
sis set in terms of operators a.i , as in the preceding 
cases. The solution was obtained by discretizing k in the 
Brillouin zone of the crystal, whose unit cell is defined 
in a coordinate system oriented at 45° with respect to 
the original axes as shown in Fig. [S] The seven distinct 
energy levels thus obtained are plotted in Fig.|S|for some 
of the main directions in fc-space. 

There are no localized states, however there is a dis- 
persionless direction, k y = 0, for excitations located on 
ribbons (or strips) which are localized in the x but ex- 
tended in the y direction. As shown in Fig. [7] the sites 
with non-zero amplitudes lie on the corners of rows of 
rectangles (formed by the double squares). Each rectan- 
gle has a pair of sites of z = 3, and a pair of sites of 
z = 4, whose amplitudes are different (corresponding to 
open (closed) circles). This state corresponds to the en- 
ergy band E w 3.85 which is, as can be seen in Fig. El 
flat in the OX direction. 

The ground state energy per site for this system is 
found within spin wave theory to be Eq/N s = —0.6475, 
close to the value obtained in spin wave theory for the 
quasiperiodic octagonal tiling, of Eq/N s — —0.646 Q. 
The quantum Monte Carlo results are —0.6602(1) for 
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FIG. 6: The seven crown lattice energy bands along some of 
the main symmetry directions. Special points in the square 
Brillouin zone are O: origin, X: zone boundary along k x axis, 
Ti: the it, n vertex. 





FIG. 8: Dependence of the local staggered magnetization 
rrii on the coordination number z for the inequivalent sites 
of the crown lattice. Quantum Monte Carlo data (open) are 
compared to linear spin wave theory results (full). 




FIG. 7: States on strips (see text). Amplitudes are different 
on the sites with open (filled) circles. 



0.36 



0.34 




0.32 



the crown lattice and = —0.6581(1) for the octagonal 
tiling 0, respectively. 

The linear spin wave theory results for the local stag- 
gered magnetization are shown in Fig. [S] for different val- 
ues of the coordination number z of the sites, and com- 
pared to the quantum Monte Carlo results. We find that 
linear spin wave theory does rather well in describing this 
inhomogeneous system. Apart from an overall overesti- 
mation of the antiferromagnetic order, both the spin wave 
theory and quantum Monte Carlo results exhibit similar 
qualitative features: for z = 3 there are two different 
values of the staggered magnetization — the sites labeled 
2 and 4 have a larger staggered magnetization than the 
site labeled 3 (more is said about this below). The case 
of z = 4 is particularly interesting, since these sites are 
found by quantum Monte Carlo to have the smallest 
value of rrii. As a result one finds a non-monotonic 
dependence of the local staggered magnetizations on z 
(comparing z — 4 to z — 5). Such a non-monotonic de- 
pendence of the local staggered magnetization on z was 
noticed in the quantum Monte Carlo |3J and numerical 
spin wave calculations for the octagonal tiling We 
will return to this phenomenon in the following section. 

For the crown lattice, the number of bonds per struc- 



FIG. 9: Dependence of the local averaged bond strength bi 
on the coordination number z for the inequivalent sites of the 
crown lattice. Quantum Monte Carlo data (open circles) are 
compared to linear spin wave theory results (full circles). 



tural unit cell is 14, but only 8 are inequivalent due to 
the mirror symmetry. From their respective strengths, 
we obtain the values of the averaged local bond strength 
shown as a function of z in Fig. EJ where the results from 
linear spin wave theory are compared to the quantum 
Monte Carlo data. We again find a difference between 
the two types of z — 3 sites. Furthermore, we observe a 
monotonic decrease of bi with the coordination number 
z. This indicates that correlations beyond the nearest- 
neighbor distance are important for the non-monotonic 
behavior of the local staggered magnetization nij seen in 
Fig.H 



IV. HEISENBERG STARS 

In this section, we present a semiclassical theory for 
Heisenberg stars, which reconciles our finding of a mono- 
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tonically decreasing rrii(z) with the intuition that increas- 
ingly highly coordinated lattices (in d > 2) should ap- 
proach a saturated order parameter. 

A Heisenberg star, studied in detail from a different 
perspective in Ref . |g, consists of a central site coupled to 
C neighbors, which are mutually disconnected (Fig. ITU)) . 
We denote the destruction operators of the Holstein- 
Primakoff bosons for the central spin, s, with a, and 
those for the outer spins, Si, with Bi, i — 1..-C- We 
then have 

C 

H = -(S 2 + Sj2 afa + B l B i + aB i + ■ ( n ) 

1=1 

Next, we carry out a unitary transformation, U , between 
the bi: hi = UijBj, with Uu = I/a/C, so that 



H 



<S 2 



(12) 



Ca f a + b\h + VC(a6i + a ] b\) + h \ b * 



i=2 



All modes hi with i > 2 correspond to precessions of outer 
spins around the central spin so that their net transverse 
component vanishes - the analogue of the local modes 
for the dice lattice mentioned above. 




FIG. 10: A Heisenberg star: a central spin is coupled to £ 
mutually un-coupled neighbors. 



The lowering of the magnetization is due to the anoma- 
lous terms linking b± with a. This part of the problem 
can be diagonalised with a standard Bogoliubov trans- 
formation matrix 



O 



( cosh0 sinh.0 
\sinh(9 cosh# 



with 



tanh(20) =2v^/(C + l) . 



(13) 



(14) 



This yields for the magnetization of the central and outer 
spins, respectively, 



m = S" - 1/(C - 1) 
M = S-1/[C(C-1)]. 



(15) 



Note the following two features. Firstly, as £ grows, the 
quantum reduction to the sublattice magnetization de- 
creases, keeping with the idea that a high exchange field 



implies a reduction in the importance of quantum fluctu- 
ations. Secondly, however, it is systematically the highly 
coordinated central site which has the lower magnetiza- 
tion: at large £, the central m is reduced by l/£, whereas 
the outer M are each only reduced by 1 /C 2 • 

The origin of this phenomenon can be traced to 
the structure of the Hamiltonian (Eq. Ill|) : Holstein- 
Primakoff bosons in a bipartite antiferromagnet are cre- 
ated only in pairs, one each on opposite ends of each 
bond. The higher a site's coordination, the more anoma- 
lous terms in the Hamiltonian act to create such boson 
pairs. As the number of bosons quantifies the reduc- 
tion in rrii, high coordination thus implies low magneti- 
zation. This is reminiscent of the order-by-disorder ef- 
fect encountered in a simple hopping problem, where one 
finds that in the ground state, particles are more likely 
to be found near highly coordinated sites. At the same 
time, a high exchange field (corresponding to large diag- 
onal terms in the Hamiltonian) makes such bosons costly 
and suppresses their overall number - which is why high- 
z lattices have a higher m than low-z ones. 

The above order-by-disorder mechanism is operative 
only when differently coordinated sites are located at op- 
posite ends of a bond. If, for the sake of argument, one 
were to consider an inhomogeneous system formed by 
patching together large pieces of uniform lattices (say, 
square and honeycomb), the local m in the interior of 
each piece would of course reflect that of the relevant par- 
ent lattice - and in this case, m(z) would increase with z. 
In other words, as placing a boson on a site involves plac- 
ing a boson on one of a neighboring sites as well, m(z) 
is of course not a purely local quantity. The next-best 
guess would be that, for sites with the same z, the mag- 
netization will be lowest for those sites whose neighbors 
have the lowest average coordination. This is in keeping 
with the result that 0.3366(5) = m 2 > m 3 = 0.3214(5) 
for the crown lattice and also with the bimodal split- 
ting of the local staggered magnetization observed for 
the five-fold coordinated sites on to quasiperiodic octag- 
onal tiling |3j,|4|. In other words, this analysis shows that 
the monotonicity in m,i(z) is a trend, not a law, in that 
details of the mode spectrum can in fact violate it. 

We close this section with two side remarks. Firstly, 
note that the Heisenberg star for £ > 1 has a non- 
vanishing total magnetization, M., which has enabled us 
to obtain a finite reduction of m due to quantum fluc- 
tuations, even though it is a 0-dimensional system. For 
Q = 1, when the the total magnetization vanishes, the 
corrections (Eq. 1151) themselves diverge. Secondly, since 
M is invariant under the action of the Hamiltonian, it 
follows that S — m = ((S — M). By the same token, 
simply dressing the Neel state on the dice lattice leads to 
S — rriA = 2(S — tub), a result which holds not only for 
linear spin waves but also in our quantum Monte Carlo 
simulations. 
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V. DISCUSSION AND CONCLUSIONS 

We have presented results for local properties of bi- 
partite inhomogeneous quantum magnets. They lead us 
to conclude that in periodic systems containing inequiv- 
alent sites, there is in general a trend towards decreased 
values of mi and hi with increasing local coordination 
number. This trend was also previously observed in a 
quasiperiodic system. The results for the crown lattice, in 
particular, show a qualitative similarity with the results 
that were earlier obtained for the self-similar quasiperi- 
odic octagonal tiling. We have provided a simple ana- 
lytical model system, the Heisenberg star, the semiclas- 
sical solution of which explains this qualitative behavior. 
Whereas it is well known that quantum fluctuations can 
lead to a range of ordering phenomena in otherwise dis- 
ordered magnets Q|, it is interesting to ask whether a 
novel disordering effect can be obtained by increasing z 
for a sublattice, as suggested by a decreasing rrii(z). How- 
ever, as we have shown above, an increase of z on one 
sublattice will generically make the ordering more robust 
everywhere, although not uniformly so for all sublattices. 
This reconciles in a simple way the observed decrease of 
rrii(z) with increasing z with the naive expectation that 
the ordered moment should grow with the coordina- 
tion of the lattice. 

SW acknowledges LPS, Orsay and CEA Saclay for hos- 
pitality during this collaboration. The numerical calcula- 
tions were performed at HRLS Stuttgart and NIC Julich. 
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